Molecular dynamics and theory for the contact values of the radial distribution 

functions of hard-disk fluid mixtures 
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We report molecular dynamics results for the contact values of the radial distribution functions 
of binary additive mixtures of hard disks. The simulation data are compared with theoretical 
predictions from expressions proposed by Jenkins and Mancini [J. Appl. Mech. 54, 27 (1987)] and 
Santos et al. [J. Chem. Phys. 117, 5785 (2002)]. Both theories agree quantitatively within a very 
small margin, which renders the former still a very useful and simple tool to work with. The latter 
(higher-order and self-consistent) theory provides a small qualitative correction for low densities and 
is superior especially in the high-density domain. 



I. INTRODUCTION 



Model systems of hard disks and hard spheres are use- 
ful for the derivation of rigorous results in statistical me- 
chanics as well as in perturbation theories of fluids. "'^ Hard 
disks and hard spheres are also relevant for the model- 
ing of mesoscopic systems such as colloidal suspensionsSi 
and granular matter^ Apart from its academic interest, 
the study of two-dimensional systems is important in the 
context of monolayer adsorption on solid surfaces?^ Re- 
cently, the equation of state of hard disks has been ex- 
perimentally measured in charge-stabilized colloidal par- 
ticles suspended in water and confined by a laser beam.^ 
While most of the studies are restricted to monodisperse 
fluids, it is obviously important to consider the poly dis- 
perse character of the system, especially in applications 
to mesoscopic matter. The equation of state, as well as 
nonequilibrium transport properties, of bidisperse sys- 
tems of inelastic hard disk s have been discussed in the 
literature, see Refs. lfiy7 y8l9yi(lllllll2yi,lll4lll,'illl(Tl and ref- 
erences therein. 



The state of an additive m-component fluid mixture 
of = X]"=i -^i hard disks is characterized by the to- 
tal number density p, the set of mole fractions x = 
{xi, X2, . . . , Xm}, with Xi = Ni/N, and the set of diam- 
eters <T = {(Ti, (T2, . . . , <Jm}- Instead of p, the area frac- 
tion V = {iT/A)p{a'^), where (cr") = J2i^i'^i is the n-th 
moment of the size distribution, can be used to charac- 
terize the density of the system. The spatial correlation 
between two disks of species i and j separated by a dis- 
tance r is measured by the radial distribution function 
(RDF) gij(r; z/, X, cr). The contact values 



Xtj [v., X, cr) = g^j {a,j ; i^, x, cr), cr^ = (ctj + CTj)/2, (1) 



of the RDFs are of special interest since they appear in 
the Enskog kinetic theory of dense fluidsii and, more 
important, they are directly related to the equation of 



state (EOS) of the fluid via the virial theoremji^ 
Z(i/,x,cr) = — ^ 

= l + 2iy^XiXj—^Xtj{'^,x,(^)- (2) 

Alternatively, the compressibility factor Z can be ex- 
pressed asiS 



Z„,(:^,x,cr) = ^XiX™('^,x, cr), 



(3) 



where Xiw{i^, x, cr) denotes the density of species i at con- 
tact with a planar hard wall, relative to the associated 
bulk density. Its expression can be obtained from that of 
Xij(z^, X, cr) by assuming the wall to be a component of 
the mixture present in zero concentration and having an 
infinite diameter>i2i2£ 



X™(!^,x, cr)= lim 



(4) 



The subscript w in has been used to emphasize that 
Eq. Q represents a route alternative to Eq. 101 to get 
the EOS of the hard-disk polydisperse fluid. Of course, 
Z = Zw in an exact description, but Z and Z^ may dif- 
fer when dealing with approximations. Thus a stringent 
consistency condition for an approximate theory of Xij is 
to yield the same EOS through both routes. 

The aim of this paper is two-fold. First, we present 
(new, accurate) molecular dynamics results for Xij in 
the case of a bidisperse hard-disk fluid mixture with a 
size ratio (Ti/a2 = 1/2. Next, those simulation data are 
compared against theoretical predictions by Jenkins and 
Mancini^ and by Santos et ali2i As will be seen, both the- 
ories agree quantitatively well with the simulation data 
but the latter is slightly superior in the high-density fluid 
regime (0.3 ^ < 0.7). The theoretical proposals for Xij 
are presented in Sec. ITU and the simulation method is out- 
lined in Sec, mil The results arc presented and discussed 
in Sec. IIVI and we close the paper in Sec. with some 
final remarks. 



II. THEORETICAL APPROXIMATIONS 

In this section, different approximations for the contact 
values of the radial distribution function are reviewed, 
first for the monodisperse and then for the polydisperse 
case. 



A. The monodisperse case 

For the sake of completeness, let us begin with the 
one-component fluid before considering the more general 
polydisperse fluid. Perhaps the most widely used ap- 
proximation for the contact value x(j^) of the RDF of 
the monodisperse hard-disk fluid is the one proposed by 
Henderson in 1975:-^ 



1 - vY 



(5) 



Despite the simplicity of this prescription, it provides 
fairly accurate values. On the other hand, Eq. (0) tends 
to overestimate the value of x{^) for high densities of the 
stable fluid phasei2i2i2^ This led Verlet and Levesque22. 
to propose the correction 



26(1 - vY 



(6) 



A more accurate prescription has recently been proposed 

by one of USi^6,27,28,29 



2^(1-^) 



(7) 



This is confirmed by Fig. ^ where simulation data 
Qf T^^j^-j^26j27j28j30 j-glative to Henderson's approximation 
X^(^'), are compared with the ratios {v) / [v) and 
X^'^ {v) / iy) . We observe that ~x^{v) behaves satisfac- 
torily up to an area fraction v k, 0.3. However, as the 
density increases and approaches the limit of stability of 
the hard-disk fluid {v ~ o.7),26,27,28,29^3i ^h^^,) overesti- 
mates the simulation data (by a few percent at v = 0.68), 
while Eq. Q presents an excellent agreement with com- 
puter simulations for densities v < 0.68. 



B. The polydisperse case 

In this subsection, the classical result for the RDF con- 
tact valued is confronted to more recent findings. The 
derivation is detailed for bulk- and wall-EOS in both 
cases, and the agreement/disagreement of the two ap- 
proaches is discussed. 
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FIG. 1: Density dependence of the ratio x/x^i where unity 
(dotted line) corresponds to Eq. 0. The circlea^£*2£i2^ and 
triangles-- are simulation data, while the dashed and solid 
lines are the theoretical predictions given by Eqs. ||SJ and 0, 
respectively. 



1. Jenkins and Mancini's approximation 

In the case of a hard-disk mixture, a useful approxima- 
tion for the contact values Xij was proposed by Jenkins 
and Mancini in 1987i& It reads 



(z^,x. 



1 



16 (1 



rzjj(x,cr), (8) 



where the parameter 



(9) 



contains the whole dependence of Xij on the size composi- 
tion through the first two moments. It is worth mention- 
ing that Eq. © was originally proposed in the context 
of inelastic disks and thus it is much better known by re- 
searchers in granular matter theory than by researchers 
in liquid theory. As a matter of fact, Eq. (jSJ has recently 
been rediscovered by liquid theorists. "^^ 

In the special case where all the disks have the same 
diameter (ai — > cr), one has — > 1 and so Eq. (jSJ re- 
duces to Henderson's approximation (jSJ for monodisperse 
disks. Thus, Jenkins and Mancini's approximation (jS)) 
represents a simple, straightforward extension of to 
the polydisperse case. As a consequence, it inherits, by 
construction, the limitations of Henderson's equation (O 
for high densities (cf. Fig.0. Moreover, Eq. ||SJ| does not 



3 



yield the same EOS through routes ^ and (PJ. Let us 
first consider the standard route Q • Taking into account 
the mathematical identities (for an arbitrary number of 
components) 



(10) 



2(^2) 

insertion of Eq. (|SJ into Eq. 10) yields 



(11) 



Z-'^^'(i/,X,(T) = 1 



l + 2iy 



\ — V (1 — z^)^ 

1-^ 



2(1 -i.) 



(12) 



where ^ = {a)'^/{a'^) was used as a convenient 

abbreviation. ^iSS To explore the alternative route 

let us take the limits indicated in Eq. on Eq. ||HJ) to 

get 



X™(^,x,^) 
where 



Zj„,(x, cr), (13) 



Zi^(x, cr) = 2a-i 



Inserting Eq. H13|l into Eq. Q and making use of 



(14) 



(15) 



one has 



Zr(^,x,^) = 5]x,x™(^,x,<t) 



= Z™(j.,x,cr)+y4 



1^/8 
(1^ 



(16) 



The inconsistency Z'"^' 7^ appears already to first 
order in v. Comparison with simulation results shows 
that is clearly superior to Z^^. 



2. Santos, Yuste, and Lopez de Ham's approximation 



The two limitations of xjf^ just mentioned, namely the 
slavery to Henderson's equation and the failure to give 
a common EOS through Eqs. Q and are remedied 



by a recent proposal made by Santos et ali2i It reads 
1 



1 - 

+ {2-u)xM- 
'1 - 1^/2 



2 - iy/2 



1 



- (1 - ^)x{^) 



Zy(x,cr) 



4- (x, cr) 



(17) 



where Zij is again given by Eq. © and the contact value 
x{v) of the monodisperse fluid can be freely chosen. Ob- 
viously, in the trivial case where all the disks have the 
same size, Zy = 1 and so ')^J^{v) — x{^)- Insertion of 
Eq. H17() into Eq. (0) yields the following simple form: 

1 



7SYH 



(i/,x,(t) = 1 



1 - 
2v 



- + Av 

J 

\-A 



Ax{y) 



1 



2(1 -z.) 

where use has been made of Eqs. (|10|l . Hll|l. and 



(18) 



E2 2 
Xi Xj (T^j Z^j 



E 



2 2 



(19) 



valid again for any number of components. Note the 
identical form of the second lines in both Eq. p2|) and 
Eq. lfTF|l . The expression for xfj^ is given by Eq. itTTIl . 
except that Zij must be replaced by Ziw When xfj^ is 
inserted into Eq. ^ and use is made of Eq. (|15|l and of 

^x,zL=4A (20) 

■i 

it turns out that the EOS (|18|1 is consistently reobtained, 

i.e. ZSYH = ZSYH. 

So far, the monodisperse quantity xi'^) remains arbi- 
trary. From that point of view, Eq. (|17ll represents a 
consistent class of approximations, with free x(j^), rather 
than a specific approximation. 

3. Some comments 

It is worth mentioning that Eqs. lO and l|17() share 
the property that, at a given packing fraction the 
whole dependence of Xij on the composition (x, cr) of 
the mixture appears through the parameter only. To 
clarify the implications of this, let us consider two mix- 
tures M and M' having the same packing fraction v but 
strongly differing in the set of mole fractions, the sizes 
of the particles, and even the number of components, 
i.e. (x, cr) ^ (x',(t'). Suppose now that there exists 
a pair ij in mixture M and another pair i' j' in mix- 
ture M' such that Zij{x,(T) ~ zti ji {s.' . cr') . Then, ac- 
cording to Eqs. © and (|17|) . the contact value of the 



RDF for the pair ij in mixture M is the same as that for 
the pair i'j' in mixture M'. This sort of "universahty" 
ansatz, which is more general than Eqs. JHl or H17|l and 
is shared by other weU-known proposals for Xij of hard- 
sphere fluids (e.g. the scaled-particle-theory, Percus- 
Yevick, and Boublfk-Grundke-Henderson-Lee-Levesque 
approximations),-^ is of course only approximate. How- 
ever, its enforcement leads to the construction of simple 
and accurate proposals for Xij with the help of only a 
few requirementsiSLS 

Interestingly enough, the EOS (|12|) and H18|l are iden- 
tical when X = is used in the latter, even though the 
contact values Xij used in the derivation are different. 
More specifically, if X = is used in Eq. ((T7|l . then 
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(1 - Zy) 



(21) 



The property Z' 



SYH 



z 



JM 



2 2 



is just a consequence of 
as follows from Eqs. 



(HU and 

If one chooses x(i^) — x^{^) then Eq. H17|) . being con- 
sistent with the condition Z = Z„,, can be expected to 
become more accurate than Eq. © , especially for highly 
asymmetric mixtures. Since x^i'^) is fairly good for 

< 0.3, as Fig. n shows, the main difference between 
Eqs. © and H17|l in that density domain lies in the func- 
tional relation on the parameter Zij : linear in the case of 
Eq. ||HJ), quadratic in the case of Eq. I|17|) . As a conse- 
quence, 



XiT / Xi 



> 1 

= 1 

< 1 if 2;. 



if < 1, 



if Zij = 1, 



> 1, 



(22) 



where we emphasize that (|22|l refers to < < 0.3. This 
is illustrated in Fig. 13 where the ratio Xi7^/Xif^ is plot- 



ted versus Zij for v — 0.25 and v — 0.5. In the former 
case, since x^ and x^^ yield practically the same value, 
the associated two functions given by Eq. (|17(l are hardly 
distinguishable. On the other hand, those functions ex- 
hibit visible small differences at v — 0.5 due to the fact 
that x^ deviates from x'"^ by about 0.5% (cf. Fig.^. 

In the spirit of the approximations H17|) and (|18|l . the 
more accurate the monodisperse function x{^) tii6 better. 
Therefore, more successful predictions for Xij and Z can 
be expected if one chooses for x{^) an expression more 
refined than Henderson's, such as Eq. iQ, especially for 
high densities. In Sec. II VI we will check all these expec- 
tations by comparing molecular dynamics results for Xij 
against Eqs. © and (|17|) . the latter being complemented 
by the monodisperse prescriptions (jSJ and Q. 



III. SIMULATION METHOD 

Since we are interested in the behavior of rigid parti- 
cles, we use an event-driven (ED) method that discretizes 
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FIG. 2: Plot of the ratio Xi/^/Xif^ against the parameter 
Zij for 1/ — 0.25 (thin lines) and v = 0.50 (thick lines). The 
dashed lines correspond to the choice = X^i'^) in Eq. 

117II . while the solid lines correspond to the choice xi'^) = 
X^'^{i^)- The dotted line represents the JM expression ©. 



the sequence of events with variable time steps for all 
particles between collisions, as adapted to the problem. 
This is different from classical molecular dynamics sim- 
ulations, where the time step is usually fixed for the nu- 
merical integration of the equations of motion. 



A. Collision Model 

The particles are assumed to be perfectly rigid and 
follow an undisturbed motion until a collision occurs as 
described below. A change in velocity can occur only 
at a collision, and due to their rigidity, the disks collide 
instantaneously. The standard interaction model for in- 
stantaneous collisions of particles with diameters (Ti and 
mass rrii is used in the following. The post-coUisional ve- 
locities (v^, v^) of two collision partners ij are given, in 
terms of the pre-collisional velocities (v^, v^), by 



2 



(23) 



mirrij / {rrii + rrij) is the reduced mass and 
Vj) • n] n is the component of the relative 
- Vj parallel to the unit vector n pointing 
along the line connecting the centers of the colliding par- 
ticles. If two particles collide, their velocities are changed 
according to Eq. (j^Sl- 



where fiij - 

v„ = [(Vj 

velocity Vj 
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B. Algorithm 



In the ED simulations, the particles follow an undis- 
turbed translational motion until an event occurs. An 
event is either the collision of two particles or the col- 
lision of one particle with a boundary of a cell (in the 
linked-cell structure) i24 The cells have no effect on the 
particle motion here; they were solely introduced to ac- 
celerate the search for future collision partners in the 
algorithm. 



Simple ED algorithms update the whole system after 
each event, a method which is straightforward, but in- 
efficient for large numbers of particles. In Ref. „35 an 
ED algorithm was introduced which updates only those 
two particles involved in the last collision. For the algo- 
rithm, a double buffering data structure is implemented, 
which contains the 'old' status and the 'new' status, each 
consisting of: time of event, positions, velocities, and 
event partners. When a collision occurs, the 'old' and 
'new' status of the participating particles are exchanged. 
Thus, the former 'new' status becomes the actual 'old' 
one, while the former 'old' status becomes the 'new' one 
and is then free for the calculation and storage of possi- 
ble future events. This seemingly complicated exchange 
of information is carried out extremely simply and fast 
by only exchanging the pointers to the 'new' and 'old' 
status respectively. Note that the 'old' status of parti- 
cle i has to be kept in memory, in order to update the 
time of the next contact, tij, of particle i with any other 
object j, if the latter, independently, changed its status 
due to a collision with yet another particle. During the 
simulation such updates may be necessary several times 
so that the predicted 'new' status has to be modified. 



The minimum of all is stored in the 'new' status 
of particle i, together with the corresponding partner j. 
Depending on the implementation, positions and veloci- 
ties after the collision can also be calculated. This would 
be a waste of computer time, since before the time i^ , 
the predicted partners i and j might be involved in sev- 
eral collisions with other particles, so that we apply a 
delayed update schemed The minimum times of event, 
i.e. the times which indicate the next event for a certain 
particle, are stored in an ordered heap tree, such that the 
next event is found at the top of the heap with a com- 
putational effort of 0(1); changing the position of one 
particle in the tree from the top to a new position needs 
O(logA^) operations. The search for possible collision 
partners is accelerated by the use of a standard linked-cell 
data structure and consumes 0(1) of numerical resources 
per particle. In total, this results in a numerical effort of 
O(A^logA^) for N particles. For a detailed description of 
the algorithm see Ref. Issl 



C. Computation of x-ij 

The results for the RDF contact values are computed 
indirectly as 



where t^^ is the average number of {ij) collisions per 
unit time and per particle of species i, A is the area of 
the system, and T = Tij = El:j'^/Nij is the temperature 
based on the kinetic energy per particle per degree of 
freedom. Note that t^^ is proportional to Nj and hence 

Xij = Xji- 

The averages are taken over a few hundred thousand 
(low density) up to several millions (high density) colli- 
sions per particle, where the first 20-30% of the simula- 
tion time is typically disregarded, so that the average is 
taken in a reasonably equilibrated state. 

IV. RESULTS AND DISCUSSION 

We have considered two hard-disk binary mixtures 
with a mole fraction of small disks Xi and a mole frac- 
tion of large disks X2 , as summarized in Table The 
diameter ratio of small to large disks is ai/a2 = 1 /2 in 
both cases. The corresponding values of the parameters 
Zij defined by Eq. © are also included in Tabled 



Mixture 


iVi 


N2 


Xl 


X2 


0-1/0-2 


Zll 


Z12 


Z22 


A 
B 


450 
7803 


126 
1998 


0.781 
0.796 


0.219 
0.204 


1/2 
1/2 


0.736 
0.747 


0.981 
0.996 


1.472 
1.494 



TABLE L Simulation parameters for mixtures A and B. 



Mixtures A and B have nearly the same composition 
(xi ~ 0.8, X2 — 0.2), but the number of particles in 
mixture B is about 17 times larger than in mixture A and 
so the statistics is better in case B. The data of Xij for 
several area fractions in the interval 0.01 < i' < 0.75 are 
given in Table Hll The values of the compressibility factor 
obtained either directly from the simulations or indirectly 
from Eq. |2Jl by inserting the Xij are also included. 

A. Contact values Xij 

In Figs. 13 E and El we have plotted the simulation 
values for the ratios Xn/xif , X12/X™, and X22/X™, 
respectively, against the area fraction i/. The data corre- 
sponding to case A are shown only for the densities not 
considered in case B, namely ly — 0.01, 0.15, 0.35, and 
0.40. The ratio Xij/Xij^ represents a "quality factor" of 
the simulation data with respect to the JM approxima- 
tion ||SJ). Figures EHSl also include the ratios xfj^^/Xij^: 
where in Eq. H17|) we have taken x(^) = X^(^) [cf- Eq. 
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TABLE 11: Molecular dynamics values of the contact values x-ij H^nd of the compressibility factor Z for the hard-disk binary 
mixtures A and B, with parameters given by table|I] The numbers in parentheses indicate the statistical error in the last digit. 



Mixture 




Xii 


Xl2 


X22 


Z 


A 


0.01 


1.013(5) 


1.015(8) 


1.02(2) 


1.0193(1) 


A 


0.05 


1.076(4) 


1.086(3) 


1.08(2) 


1.1026(3) 


A 


0.10 


1.162(7) 


1.181(5) 


1.20(2) 


1.2229(7) 


A 


0.15 


1.264(5) 


1.289(6) 


1.34(1) 


1.3656(9) 


A 


0.20 


1.383(4) 


1.426(6) 


1.48(2) 


1.537(2) 


A 


0.25 


1.518(2) 


1.579(6) 


1.66(2) 


1.742(2) 


A 


0.30 


1.689(3) 


1.768(8) 


1.88(2) 


1.997(2) 


A 


0.35 


1.888(5) 


1.998(7) 


2.15(3) 


2.310(3) 


A 


0.40 


2.132(6) 


2.28(2) 


2.51(5) 


2.708(7) 


A 


0.45 


2.441(8) 


2.63(3) 


2.98(7) 


3.22(1) 


A 


0.50 


2.828(5) 


3.09(2) 


3.52(4) 


3.89(1) 


A 


0.55 


3.346(6) 


3.68(1) 


4.35(6) 


4.79(1) 


A 


0.60 


4.03(1) 


4.51(2) 


5.41(7) 


6.05(2) 


A 


0.65 


5.01(1) 


5.65(3) 


7.03(9) 


7.87(3) 


A 


0.70 


6.47(2) 


7.36(7) 


9.6(2) 


10.67(5) 


A 


0.75 


8.77(1) 


10.14(9) 


13.6(5) 


15.3(2) 


B 


0.05 


1.076(1) 


1.084(3) 


1.09(1) 


1.1028(2) 


B 


0.10 


1.164(1) 


1.180(3) 


1.21(1) 


1.2234(3) 


B 


0.1482 


1.262(1) 


1.288(1) 


1.33(1) 


1.3609(5) 


B 


0.20 


1.385(2) 


1.423(2) 


1.50(1) 


1.5376(7) 


B 


0.25 


1.525(1) 


1.582(3) 


1.68(2) 


1.745(1) 


B 


0.30 


1.690(1) 


1.768(2) 


1.90(1) 


1.996(2) 


B 


0.3455 


1.873(1) 


1.978(5) 


2.15(2) 


2.282(2) 


B 


0.3867 


2.070(3) 


2.196(6) 


2.45(2) 


2.595(3) 


B 


0.45 


2.450(2) 


2.638(7) 


3.00(2) 


3.221(5) 


B 


0.50 


2.843(2) 


3.096(4) 


3.59(2) 


3.891(2) 


B 


0.55 


3.357(3) 


3.705(7) 


4.37(3) 


4.794(2) 


B 


0.60 


4 058f2^ 


4 525^6^ 


5.49(2) 


6 055^3^ 


B 


0.65 


5.036(7) 


5.708(9) 


7.06(5) 


7.885(3) 


B 


0.66 


5.281(6) 


6.00(2) 


7.46(5) 


8.348(4) 


B 


0.67 


5.547(4) 


6.32(1) 


7.91(4) 


8.855(1) 


B 


0.68 


5.830(5) 


6.68(1) 


8.36(4) 


9.409(1) 


B 


0.69 


6.145(8) 


7.06(2) 


8.89(5) 


10.016(2) 


B 


0.70 


6.498(9) 


7.45(3) 


9.55(7) 


10.684(3) 


B 


0.75 


8.82(2) 


10.34(4) 


13.4(1) 


15.30(1) 



(O] and x(i/) = X^'^i'^) [cf- Eq. {T))] for the monodisperse 
fluid. 



1. Low densities 

We observe that xfj^^/Xij^ with both prescriptions 
xi'^) = X^i'^) E^nd xi'^) — X^'^i'^) are practically indistin- 
guishable up to « 0.3. This is consistent with the fact 
that in that domain of low and moderate densities the 
correction to Henderson's EOS is irrelevant, as shown 
in Fig. n On the other hand, some limitations of Eq. © 
are already apparent in the range < < 0.3: the JM 
approximation (slightly) underestimates the small-small 
contact value (x'™ < Xii)i while it overestimates the 
large-large contact value (x22^ > X22)- These two effects, 
which are not linked to the use of x^('^)) *re reasonably 
well captured by Eq. (|17|l . as expected from H22|l . In the 
case of the cross contact value, we have Xi2 — X™ for 
< < 0.3, in agreement with the fact that zi2 ~ 1 in 



our mixtures A and B. 



2. High densities 

In the high-density domain v > 0.3, the simulation 
data clearly deviate from both Eq. (jS)) and Eq. H17|) when 
the latter is combined with xi^^) = X^i'^)- Both the- 
ories tend to overestimate the contact values, what is 
essentially a trait inherited from Henderson's EOS. On 
the other hand, a much better agreement is obtained 
when Eq. (|17|l is combined with Eq. Q for the monodis- 
perse fluid. The remaining deviations of the latter the- 
ory from the simulation values are a reflection of the ap- 
proximate character of Eq. (|17|l rather than that of Eq. 
{Tj), in view of Fig. ^ Part of the deviations for the 
highest densities may be due to the proximity to crys- 
tallization. In the monodisperse case, it is known that 
the hard-disk fluid undergoes a freezing transition (pos- 
sibly mediated by a hexatic phased) at an area fraction 
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FIG. 3: Density dependence of the ratio Xn/Xii* for a hard- 
disk binary mixture with parameters as given in Tabled The 
filled squares and open circles represent our simulation data A 
and B, respectively. The lines represent the ratio Xvl^ IX™ 
for case B with x(^) = 'X^i^) (dashed line) and with x(^) — 
x'^^(i^) (solid line); the dotted line represents Eq. The 
theoretical curves for case A are practically indistinguishable 
from those for case B and so they are not plotted. 




FIG. 5: Same as in Fig. except for X22/X2 



V ~ 0.7i ^?i^'i'i^?i^?i?ii?? In any case, polydispersity tends 
to increase the freezing transition density. For mixture A, 
our simulations indicate that the transition takes place 
between u = 0.75 and v = 0.8. For more details on very 
high densities in the mono- and bi-disperse situations see 
Refs.l26landl27l respectively. 



1.01 




0.0 0.2 0.4 0.6 0.8 



V 

FIG. 4: Same as in Fig.El except for xri/xlf- 



B. Equation of State 

Although in this paper we have been mainly concerned 
with the contact values of the RDFs, it is worth consid- 
ering the compressibility factor Z = p/ pUbT. The simu- 
lation values of the ratio Z/Z'^^, where Z'^^ is given by 
Eq. (inji are plotted in Fig. We observe that the JM 
EOS is fairly accurate for v < 0.4, even though the indi- 
vidual values Xij^ ^^^e not that good in the same density 
range. This is mainly due to a fortunate "cancellation of 
errors" (x™ < Xiij while > X22)- As a matter of 
fact, as already mentioned in Sec. m the recipe H18|l in 
combination with x(^) = X^i'^) becomes identical with 
Z"^^. Nevertheless, the JM approximation again overes- 
timates the simulation data for higher densities {ly > 0.4). 
When x(i') = X^'*(^) is used, Eq. lfTH|l becomes quite rea- 
sonable, although it slightly underestimates the mixture 
pressure of the fluid at the highest densities. 

Interestingly, a recently proposed empirical correction 
to the EOS 

Z = 1 + (1 - aiy^) (Z™ - 1) , a = 0.1, (25) 

see Eq. (20) in Ref. |29|, works also pretty well for the 
parameter set used here, however, without theoretical 
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FIG. 6: Density dependence of the ratio Z/Z''^ for a hard- 
disk binary mixture with parameters as given in Table |I] The 
filled squares and open circles represent our simulation data A 
and B, respectively. The lines represent the ratio ^syh^^.tm 
for case B with xi'^) ~ X^i'^) (dotted line - unity) and xi'^) ~ 
X^'^i^) (solid line). The dashed line is the empirical relation 
(Unj proposed in Ref. ^ 



foundation. The excellent performance of Eq. (|25|) for 
xi ~ 0.8 and iTi/ct2 = \ does not necessarily extend to 
other compositions. In fact, Eq. H25I) is not as good as 
Eq. 10 in the monodisperse case. 

V. CONCLUSION 

In this paper we have presented molecular dynamics 
results for the contact values Xij = 9ij {<^ij ) of the radial 
distribution functions (RDFs) gij{r) of binary mixtures 
of additive hard disks. As a representative case we have 
fixed mole fractions xi ~ 0.8 and X2 — 0.2, and used a 
diameter ratio ai/a2 = \- A set of numerical values of 
the area fraction v have been considered, covering the 
dilute, the intermediate, and the dense regime. 

The simulation results have been used to assess the re- 
liability of theoretical expressions previously proposed in 
the literature. Until recently, practically the only pro- 
posal was the one by Jenkins and Mancini,^ which is 



expressed by Eq. ((SJ. This approximation succeeds in 
capturing the main trends in the intricate dependence 
of Xij on the area fraction and the composition of the 
mixture, even at a quantitative level. 

On the other hand, our simulation data expose some 
(small) limitations of xlf^' already for low and moderate 
densities {v < 0.3) the JM recipe tends to underesti- 
mate the small-small correlation value and overestimate 
the large-large value; for higher densities {v > 0.3), Xif^ 
overestimates the correct values, an effect that can be 
traced back to the fact that xfj^ is strongly tied to Hen- 
derson's EOSiSa 

These two shortcomings are widely corrected by a re- 

Eq. in|). While 



cent proposal made by Santos et al.pi Eq, 



Xij^ is a linear function of the parameter Zy defined by 
Eq. (O, Xi/^ is a quadratic function. This higher or- 
der approach allows xfj^^ to satisfy an extra consistency 
condition in the limit of highly asymmetric mixtures^ 
Moreover, xf/^ can be used in conjunction with any de- 
sired expression for the contact value of the monodis- 
perse fluid. When instead of Henderson's expression the 
one recently proposed by one of us^SiSi is employed, Xi/^ 
exhibits a reasonable agreement with the simulation data 
for 1/ > 0.3. 

In spite of this, it can be observed that xf]^^ tends 
to underestimate the simulation data for very high den- 
sities (v > 0.6), so that an even better approximation 
is needed in that extreme, high density fluid regime. 
From this point of view, we hope that our simulation 
data will be helpful to test the accuracy of other future 
theoretical proposals that have been or will be made. 
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